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S YSTBM AND METHOD FOR OPTIMIZING DRUG THERAPY FOR THE 
TREATMENT OP DISEASES 

This application claims .priority benefit of European Patent Application No. 
5 00/203200.1, filed on September 15, 2000, and U.S. Provisional Application No. 
60/279674 March 30, 2001, the contents of which arc expressly incorporated by 
reference herein. 

Field of the Invention 
10 The present invention generally relates to the field of drug therapy, disease 

* managem ent therapy monitoring and phannaoogenoniicsL In rae embodiment, the 
present invention idatea ix> systems 

therapy for a patient in connection with the treatment of a disease. The present 
invention also provides an approach towards therapy design based on the integration of 
IS bio-analysis, pharmacological modelling and resistance testing. 

Bactoasw^.tf % foreafrm, 

Infectious agents including tuberculosis bacillus, human immunodx£dency virus (HIV), 
and cell proliferative dis codex* have proven difficult to treat once aflecting an 

20 individual. Efficacy of antrretro viral therapy is generally measured by a drop in viral 
load (concentratioii of viral RNA copies in the blood plasma), while aofoBtroviral 
therapy failure is generally reflected by an increase in viral load and/or the 
development of resistance to therapy likewise, anti-cancer drog treatments and 
therapies (Le, f chemotherapy, gene therapy, radiation, etc) have proven effective 

25 against many malignancies and forms of cancer. However, many patients experience 
treatment failure, or reduced efficacy oyer time with many anti-cancer drugs and 
therapies. Such treatment failure may be due to a variety of causes, such as 
development of resistance to the particular drug via mutation or other process, 
progression of disease requiring an altered dosage regimen, patient noncomp liance, 

30 sub-optimal pharmacokinetics, toxicity to a drug etc. 

Intermittent blood level m oni t orin g of drugs has been described in die literature as 
'therapeutic drug monitoring. 1 ' True theraporic drug jnom 

accurate, would require constant, quantitative drug monitoring of blood concentrations 
35 in each individual patient for each administered drag. However, besides being 
prohibitively invasive and time consuming, such an approach suffers firm various 
other practical shortcomings. Since such actual, constant blood level -monitoring of all 
administered drugs is nearly impossible, some interval between samplings is required; 

CONFIRMATION COPY 
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different drn^ maybe adntimstercd at different times post^dministr^tion, leading to 
irregular sampled drug disposition carves. ; 

Treatment ffiiccffffl for ^'™« t inftlh^^E cancer, infectious diseases and viral 
5 illnesses, is correlated with the use of optimal drug dosages, for both single drugs and 
for drugs in combination. Optimal dosages guarantee that the plasma drug 
concentration^) remain well above the mrnitTmm effective concentrations (MECs) of 
aU the arhmmstei^ of a particular drag 

in a particular patient, the to wear flic disease sensitivity is to that particular drag, 

10 resul ting in lower likelihood of effective treatment Hie j^bability of treatment 

success dep ends on the feet that the MEC is drug-specific, and that for the same drug ,. 
the MEC also varies across the patient population. Abo, different drugs are mora 
effective in some patients than in other patients due to inter-individual differences in 
phar macokineti cs. Individual patient characteristics also effect dosages, Le^ 

15 characteristics such as body size, gender, age, physical and pathophysiologic states, 
genetics, environment, and conc u rr e nt therapies. Therefore current day therapeutic 
monitorning services based on the sole detennmation of the concentration of a drug in a 
sample of a patient may have limited value. 

20 Previous research has attempted to navigate effective dosages of drugs to c ha l l eng e 
rapidly changing etiologfc agents. While the broad approach of population 
pharmacokinetics (loosely defined as the change in time of thq concentration or nature 
of therapeutic agent(s) in groups of patients having similar characteristics) is a 
technique of long standing (see t_M_ Ladder^ J. CHn. Pharmacol, 28:1059-1062 

25 (1988)), it fails to take into account a large amount of inter*, and even inter, patient 
variability, ultimately contributing to therapy failure. This is in part completed by the 
development of Bayesdan parameter estimation in conjunction with population 
phannacokmetics (Thomson & Whitinfr Clin. Pbatmacokmet; 1 992, 22(6% 447-467). 
The combination of these parameters provides an approach to determine patient 

30 specific pharmacokinetic variables. 

Another difficulty in the field of drug therapy is the development of drug resistance, 
which further stresses the need for individualized therapy. Fox example, corrtinnous 
high level in vto replication of retroviruses, particularly HIV, and the intrinsic error 
3 5 rate of the reverse transcriptase enzyme are-major driving forces behind the generation 
of drag resistant vims variants. When sub-optrmal drug dosages are applied as a 
pressure to tins divergent and rapidly replicating virus population, variants with fee 
appropriate mutations in their genome will escape drug inhibition and outgrow the 
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wild-type, drug-suscej^ 

are faced with ever narrowing therapeutic options. HIV drug resistance is an ever 
increasing problem, with an estimated 10 to 20% of patients in developed countries 
frrilwig to respond to drag therapy in the first year of treatment and dev doping 
5 resistance to at least one drug. 

Likewise, malignant ceDs, such as tumor cells, are subject to similar selection pressure 
by sub-optimal dosage therapy. Mutations accumulate ova- time, resulting in malig- 
nancies recalcitrant to drug therapy. One example of a specific mutational target is the 

10 tumor suppressor gene p53. The tumor suppressor gene p53, located on chromosome 
17, is a key component, of the body's antitumor defense (Soossi, T\; Aim. -W.H Acad 
Sci„ 910:121-139 (2000); North, S, & Hainaut P.; Pa&oL Biol 48:255-270 (2000); 
Somasundaram, XL; Front Bios*. 5: D424-437 (2000); Tokino, T. & Natamm, Y.r 
Crit Rev. Oncol HemttoL 33:1-6 (2000)). The p53 gene normally responds to UNA 

15 damage that might otherwise lead to cancer by arresting cell growth, initinrrng DNA 
repair, or sending cells into apoptosis (programmed cell death)* When a p53 gene is 
mutated, however, thep53 gene, and the cells expressing it, become an etiological 
agents for cancer. Not only are tumor suppressor effects lost, but uncontrolled cell 
growth is promoted, leading to increased cett d^ 

20 increases in mutation rate, and thus further cancers* As a result, an imirvidual patient's 
resistance to available treatments (e.g., cancer treatment, antiviral therapy) also must be 
taken into account when determining an effective therapy regimen, 

r>rug resistance, or thcra^ 

25 genotypic testing, or by a combination theroo£ Drug resistance^ or therapy resistance, 
is generally determined by two main methods, namely phenotypic testing and 
genotypic testing, or by a combination thereof Phenotypic testing directly measures 
lb* »<*im \\ thmqy wig*™™* f>fa pnriftnt*g Tnntignimt nr infected eells to a particular 
therapy or therapies (generating, for example, a concentration of that drug which results 

30 in a 50% inhibition of virus gmwth, Le-, the IC50). The phenotypic testing measores 
tiie ability of a vims, for example, to grow in the presence of various drugs. Genotypic 
resistance testing (sometimes called genotyping) examines the genetic material offlie 
cell or virus to detect the presence of specific genetic mutations or patterns of 
mutations in the gene or genes of interest that confer resistance to a certain therapy ear 

35 therapies. Genotypmg can be more rapid and less expensive than phenotyping, but may 
be more difficult to accurately interpret, due to the hundreds of mutations involved, for 
example, in HTV or p53 oncogenesis. 
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Allhcni^ j4cm>typic testing 
assessmexjt of thei^y resistant 

longer and may generally be more expensive than genotypic testing Compared with 
phenotypic testing, genotypic testing has advantages, including the relative simplicity, 
5 low cost, and the speed with which the test can be performed. Currently, genotypic 
interpretation, has predominantly been applied to determining r esis t an ce of a vims, e.g., 
HIV, or mutations in a viral stain to a therapy. In a farther development this analysis 
can be performed using the approach of virtual phenotyping (eg. VirtaalPbenotype, 
PCT/EPO 1/04445), wherein the sequence of an etiologic agent is compared to 
10 sequences present in a database. Hie corresponding phenotype can be calculated based 
on the phenotypic data, of the similar sequences. 

In addition, a therapy can be less effective or ineffective in an individual because of 
allelic variations at genes important for the action of a drug. This allelic variation can 
15 mfun ynrintion at thr drug targ# hrrt alftfl ft* C^*" ^"flneficiiig the drug pihairnacntintjcs 
and pharmacodynamics. Genes which metabolize the drug or receptors influencing the • 
distribution of said drug* 

Therefore, because of the importance of maintaming an effective MEC in order to 
20 avoid the development of disease resistance, and the need to consider an individual 

patient's resistance to known therapies in the calculation of optimal dosage of a therapy 
regime for that patient, there exists a strong need in the art Jbr a single therapeutic 
procedure to aid doctors with optimizing treatment of these diseases. There also exists 
in the art a strong need for individualized therapies and qptinrizalian of these therapies 
25 for individual patients. This need is particularly strong in view of the plasticity the 

drug response of diseases such as virus infections and malignancies* This optimization 
should be adaptable to singl e dro ffl as well as to combinations of drugs sod treatment 
and should provide a mod^l with inputs for actual individual patient data as 
well as overall population data from patients (such as from clinical trials), in order to 
30 assess tea* all known therapies whether plasma levels remain above the MEC 
throughout therapy on a patient by patient basis. 

In the art individual methods are disclosed to determine re si stance (e.g. 
Antrvirogram®), to determine &c concentration of agents in a biological sample (eg, 
35 high pressure liquid chromatography, mass spectrometry) and to model the 

pharmaco kinetics of drugs administered to individuals. Karlsson MO, Sheiner LB., 
J Fharmacokinet Binpharm 1993,21:735-750; Mandema JW, Verotta D, Sheiner LB-, 
J Pharmacokinet Biopharm 1992^0^11-528; Thomson AH. Whiting B., CHn 
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Fharmacokinet 1992^2:447-467; Wakefield J, Riuaim-Poan A., Stat Med 
1495,14:971486; Roaner OL, Mufler P., J Phannacokmet Biophaxm 1997,25:209-233; 
Bennett JE, Wakefield JCL, J FhannacokixictBicFpbflnn 1996,24:403-432. Though 
these methods provide infonnation on either variable, the individual parameters allow 
5 limited imaging patient treatment For instance, the drug level in the circulation will 
not provide evidence regarding the occurrence of resistance. The need for additional 
data apart from either drug monitoring of RN A testing in the follow-up of HIV therapy 
was described by Duranl and coworkers (AIDS, 2000, 14, 1333-1339)* This group 
linked the RNA levels to fee plasma drug concentrations. However, titi$ group did 

10 neither use population baaed modeling, nor pbenotypic data, nor the combination 
thereof to evaluate drug effectiveness. Therefore, in aider to design a therapy for 
diseases such as cancer and retroviral infections, disease states in which resistance 
displays a critical role, an integrated approach combining resistance testing, bio- 
analysis and phannacolrinetic modelling is needed to provide a patient specific therapy 

15 management This integrated approach is the subject of the instant invention. 

The present invention adds to the art a combination of a bio-analytical method with 
population based modeling to determine a patient specific measure of therapy exposure, 
and a resistance determination. The combination of the resistance and patient specific 
20 pharmacokinetic parameters provides a single measure to manage therapy. Una single 
variable provides the treating physician with a measure of therapy efficacy and to draw 
conclusions on an patient specific basis fin: cither ding dosages and resistance patterns, • 



25 The present invention relates to methods of measuring the efficacy of at least one 
therapeutic agent comprising a combination of a patient's exposure to a therapy and 
resistance data. For example, in one embodiment, the invention relates to a method of 
measuring the efficacy of at least one therapeutic agent comprising: detcmiininga 
pharmacologic exposure either using a measured or predicted population 

30 pharmacokinetic model for said at least one therapeutic agent; determining resistance of 

* 

an etiologic agent towards said at least one therapeutic agent; determining the 
inhibitory quotient for said at least one therapeutic agent based on said pharmacologic 
exposure and said resistance; and using said inhibitory quotient to determine efficacy of 
said at least one therapeutic agent In one embodiment, the methods of the invention 
35 further comprise the use of a bioanalytical method to obtain an actual concentration of 
at least one therapeutic agent ma patient The inhibitory quotient may also, for 
example, be normalized to one embodiment, die population pharmacokinetic model for 
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use many of the embochments of the invention may be moptmrized population, 
phannacddnctic modeL 

In one embodiment, fee inhibitory quotient used in practicing any aspect of fee 
5 invention w«y, for example, be determined by a method comprising 

a) obtaining' an actual concentration of at least one therapeutic agent in a 
patient at a givai time using a bionalytical method; 

b) calculating a theoretical concentration of said at least one therapeutic agent 
in said patient at said time using a first population phannacoldnetic modd; 

10 c) obtaining a difference by conujaring the theoretical concentration of said at 

least one toenttjeutic agent wife fee actual coi^tn^onofsdd at least one theraiwtfc 
agent in a patient; 

d) miinnnang the ■Wfemnce bv changing at least one parameter in fee first 
population pharmacokinetic model in order to generate an optimized population 

15 pharmacokinetic model; 

e) obtaining resistance data from said patient; 

f) detfsnuningthehfen^ 

based on said optimized population phamaacokinetic modd and said resistance. The 

mefeodinay further comprise 

20 

The inhibitory quotient, may, for example, be used to optimize at least one of a 
• merapentic agent regime, mctodmg, but not limited to fee choice of therapeutic agent, 
including combinations of therapeutic agents, and fee dosage of a feerapentic agent 

25 The invention encompasses any method or methods of generating resistance data, 
whether based on genotype, phenotype, or some combination thereof. 

The present invention also relates to inetbndsofoptnniziagatleatf ^ 

agent regime for at least one patient comprising a combination of a pharmacokinetic 

30 modd and nsistancetlata. FOTexarimle,mone i^bodiment, thcinveationrdates to a 
method ofcpthmzmg at least one feerapentic agent regmie comprising: dttommmg a 
phmmacologic exposure using a population pharmacokinetic model tor at least one 
therapeutic agent; deternmung resistance of an etiologic agent towards said at least one 
therapeutic agent; determining fee hfejbitory quotient for said at least one therapeolic 

35 agent based on said pharmacologic exposure and said resistance, and using said 
iimibitory quotient to optimize said at least one therapeutic agent regime, la one 
embodiment, fee methods of the iimaiuan former comprise the use of a bioanalytical 
method to obtain an actual concentration of at least one fhenmeulic agent in a patient 
The inhibitory quotient may also, for example, be norraahzed. 
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The present invention also relaies to methods for obtaining a dosage regime for at least 
one therapeutic agent for at least cone patient comprising a combination of a 
phannacotenetic model and resistance data. For example, in one embodiment, the 
5 invention relates to a method for determining a dosage regime for at least one 
therapeutic agent comprisirjg: determinipg a pharmacologic exposure using a 
population pharmacokinetic model for at least one therapeutic agent; determining 
resistance of an etiologic agent towards said at least one therapeutic agent; determining 
the inhibitory quotient for said at least one therapeutic agent based on said 
10 pharmacologic exposure and add resistance, and using said inhibitory quotient to 

determine a dosage regime for at least one therapeutic agent In one embodiment, the 
methods of the invention further comprise the use of a bioanalytical method to obtain 
an actual concentration of at least one therapeutic agent in a patient The inhibitory 
quotient may also, for example, be normalized. 

15 

The present invention abo relates to methods for providing advice to a physician 
regarding at least one therapeutic agent for at least one patient enrnprisnng a 
combination of a pharmacokinetic model and resistance data. For example, in one 
embodiment; the invention relates to a method for providing advice to a physician 

20 regarding at least one therapeutic agent for at least one patient comprising: determining 
a pharmacologic exposure using a population phrnmacokinetic model for said at least 
one therapeutic agent; determining resistance of an etiologic agent towards said at least 
one therapeutic agent; determining the inhibitory quotient for said at least one 
therapeutic agent based on said pharmacologic exposure and said resistance, and using 

25 said mhibitory quotient to provide advice to a physician regarding at least one 
therapeutic agent for at least one patient In one embodiment, foe methods of the 
invention further comprise foe use of a bioanalytical method to obtain an actual 
concentration of at least one therapeutic agent in a patient The mhibitory quotient may 
also, for cacamplcy be normalized 

30 

The present invention also relates to methods for providing a report regarding at least 
one therapeutic agent For example, in one m^nt^ foe invention, relates to a 

method for providing a report 00 ootp rismgp jxmnng a pharmacologic exposure using 
a population p^ |f,TT| r iaflf |i VfiiBtic mo de l ifor said a t least one therapeutic agent^ dens nnpxng 
35 resistance of an etiologic agent towards said at least one therapeutic agent; detemdning 
foe inhibitory quotient for said at least one therapeutic agent based on said 
pharmacologic exposure and said assistance, and providing a report comprising at least 
one entry chosen from foe inhibitory -quotient and information derived from the 
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inhibitory quotient Inowanlttdimentth^ 

the use of a taoanarytical method to obtain an actual concentration of at least .one 
therapeutic agent in .a -patient 3he inhibitory quotient may also, for example, be 
normalised- The invention also mcludes, for example, a report ccmprismg the 
5 inhibiioiy quotient 

In another ensbodiment, the invention relates to a computer system comprising at least 
one database comprising at least one inhibitory quotient for at least one patient The at 
least one inhibitory quotient may, for example, be a normalized inhibitory quotient 

10 , ' 

In another embodiment, the invention relates to a method of identifying at least one 
therapeutic agent effective against at least one etiological agent comprising: 
determining a pharmacologic exposure using a population pharmacokinetic model far 
said at 1 east one therapeutic agent; detenmnmg resistance of said etiologic agent 

15 towards said at least one therapeutic agent; detenzmnng the inhibitory quotient for said 
at least one therapeutic agent based on said pharmacologic exposure and said 
resistance, and using said inhibitory quotient to identify at least one therapeutic agent 
effective Ag^n^ at least one etiological agent In one embodiment, the methods of* the 
invention further comprise uieuse of a bioanalytical method to obtirin art actual 

20 concentration of at least one therapeutic agent in a patient The inhibitory quotient may 
also, for example, be npmafced. 

In a further embodiment the invention relates to a method of identifying toxic effects 
of at least one therapeutic agent cwinraring; detenxumng a pharmacologic exposure 

25 using a population pharmacokinetic model for said at least one therapeutic agent; 

determining resistance of an etiologic agent towards said at least one therapeutic agent; 
determining the rnhibitory quotient said at least one therapeutic agent based on said 
pharmacologic exposure and said resistance, and using said inhibiioiy quotient to 
identify toxic effects of the least one uierapcuhe agent In one einbodixnent, the 

30 methods of the invents 

an actual concentrati Tne inhibitory 

quotient may also, far example, be normalized. 

The invention further relates to systems, connjuterpiogn^ 
35 methods, servor aide and cUeru^ side systems and methods fur ge ner a tin g , providing, and 
transmitting optimal dosage regimens for an individual patient 

Both the foregoing general description and the following detailed description are 
exemplary and are intended to provide further explanation of the invention as claimed 
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Brief Q^criDtiaa of the Drffiw^ 

The accompanytog drawings provide a further undfsfitrodiiig of the invention and are 
incorporated in and constitute a part of this specification. The dnrwm^, together with 
the description, illustrate various anbodimetxts of the mveatian, In the drawings: 

Figure 1 is an exeo^laiy gr^ oi*8ie »^ 
time; Figure 2 is an exemplary flow chart far optimizing a therapy, in accordance 
with the methods of tie invention; 

figure 3 is an exemplary representation of a system environmoit in which the 
features and methods of the invention may be implemttted; 

Figure 4 is thoiriarionship between ampreaavirNIQ and change in viral load at 
wedc24. Circles are actual values and the lw is fitted value from 
Emax model. 



10 



20 



15 Detailed Description of f hc favention , 

These terms as used hodn are defined as follows: 

"Bioanalydcal method" or bioanalytical testing means any analytical technique 
known in the art to determine the presence and/or the amount or concentration of a 
trjoapymapati^samnle. Techniques include, but are not Iimhxd to, hi^i 
performance Hquid chromatography, mass spectrometry, LC-MS, iadioimmnnoassay, 
enzyme linked immunosorbent assay, and ou^^ 

A '"biological sample" is any material obtained from a patient which ffiVt itfit f f^ff 
etioiogicaJ agent ameaabte to therapy resistance testing. Some examples are saliva, 
semen, breast mflk, blood, plasma, feces, urines tissue samples, cells in cdl cnrfaire, 
25 cells which may be farther Stared, etc For example, in a patient infected with HTV, 
any biological sample containing virus may be used. For a cancer patient, a sample 
would include all of the above, and tumors, biopsy tis^ 
of p53 could be detenninecL 

■Clinical data* may include previously recorded patient data, innhirfing 
genotypic variation or patterns with specific therapy sensitivities, data from phcaotype- 
genotype relational databases, 50% uniilritory concentrations and mirnnmm effective 
concentrations of various therapies, known drug-drug interactions, indications, or 
coiihTuudicatiGn^ 

obtained from public databases or journals, or forwarded by rcseaxdias in the field 
35 A "copgromcalifla channel* 1 is any channel which allows cognntunieaiion 

between different people, computes, or locations, i e^ telephone lines, wireless 
networks, compute networks, public networks (such as the rntemet), private networks 
(such as an intranet), sateffite-based networks, manual entry of data into a common • 



30 
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database, etc* Una communication rJmmirf may bo digital or analog! re a l tin** or 
delayed, and oneway or two way, or any combination or combinations ihereofbetwecn 
the ^ifif ^j cat entities. 

The teem "doctors" or "physicians" is understood to include any professional 
5 person authorized or trained to treat or take patient data and/or samples. • Such persons 
include but are not limited to clinicians, health care workers, nurses, technicians, etc 

"Dosage" includes the size, frequency, formulation, comedi cation,, and number 
of doses of at least one therapy to be given to a patient This al s o includes newly 
prescribed therapies and/or therapies, both singly and in combination and is irrespective 
10 of the way of adrnhrMration. 

■Resistance" or Therapy resistance" includes any condition by winch the cells, 
etiologic agent or patient respond or adapt to a therapy. 

An "etiological agent* is a disease producing agent Examples of rapidly 
mutating etiological agents are viruses such as retroviruses, and cancer causing genes 
15 or gene mutations such as those found inp53 and other oncogenes Other agents 
include bacteria, viruses, prions, algae, fungi, and protozoa. 

"Genotypic resistance" comprises changes in the genome of a ceE, virus, or 
diseased ceU'associated with the resistance to a therapeutic agent or therapy, A 
diseased cell includes, but is not limited to, cells ™far*«i hy j» ^frns, or a ba rtrt frrn t 
20 and cells with an altered phenotype hy prolifera t i on, inflammation t ^enrratfoi 
"Genotypic testing" analyzes part or all of a genetic sequence. Tins method 
may include fill] or partial genomic sequencing by all known means, and may be 
correlated with phenotype. One such method is the Vntualphenotype® 
(PCT/EP01/04445). 

25 "EFT is the human immunodeficiency virus, which is a retrovirus and of winch 

different species are currently known. A rotavirus includes is any RNA virus that 
utilizes reverse transcriptase during its life cycle. 

50% inhibitory concentration, or ICso, is the auoount of a substance required to 
inhibit growth in 50% of cells or organisms in vitro. 
30 nnhflritay quotient", IQ, is a ratio of a measure of therapy exposure and a 

measure of viral susce^ For example, IQ is the divided by 

the ICso for a particular therapy » 

A "patient" is any organism, particularly a human ex other Tnnrmniiij suffering 
from a disease, m need of testing or 

35 screening for a disease, A patient nifhrffp **ry ipftmrn ^ n ttludiag -farm animals or 
pets, and includes banana of any age or state of development 

"Patient date* includes, but is not limited to, age, gender, weight height, 
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allergies, other therapies, physical condition, diseases states), medications 
currently bong taken, disease status or progression, etc. 

"Population phmmucokmt^ic modd" or a pharjnacokinetic model predicts an 
individual plasma concentration of a thoapeutic agent using a set of mathematical 
5 equation, An "optimized" population pharmacokinetic model is a modd that has been 
adjusted to minimize the difference between at least one data point in fee model and at 
least one actual measurement from a patient The pharmacokinetic model which 
describes the drug's behaviour in an organism can be chosen out of variety of models 
known to the person skilled in the ait, including, but not limited to, models based on 

10 one compart ment, two or more compartments, and using either zero order, first order 
second order or higher order kinetics. The model may be a predicted model, wherein 
the model is chosen based on data known in the art for a therapy. Alternatively, the 
model may be measured by analyzing patient sample and determining the 
pharmacokinetic model thereon (measured model). 

1 5 For example, based on literature data and/or drug concentration determinations 

in patient indications for a model may be provided. A model may allow one to predict 
or estimate parameters required, ag. Ctrough, Patient parameters may also be in cl u d e d 
in the modd, eg. age, gender, weight, body mass index (Bayes approach). In one 
embodiment, this combination of data and mathematic equations allows the prediction 

20 of parameters induding the dosage regimen needed to obtain a certain drug 
concentration. 

Pharmacologic exposure^ is the extent to which a patient is exposed to a 
therapy. A measure of exposure is, e.g. Ctrtmgh and area under the curve (AUC). 

"Pheno typic resistance" comprises fold-resistance compared to a reference of a 
25 cell, virus, or viraHy infected cell to a tested theraponic agent or therapy, specifically, . 
traits that can be observed Thenotypic testing 11 is a testing method that obtains tins 
trait o£ for example, a cell line or virus. One such method is the high feroughpul vinil 
screen Antivixogram<3> (Virco, Belgium; WO97/27480; US 6^221^78). 

"Etiological agent" includes any agent which causes disease in a patient Some 
30 examples include, but are not limited to viruses, particularity HIV, bacteria, and 
mutations associated with malignancies, such as p53„ 

A 'therapeutic agent' 1 is a dni& pharmaceutical, antiviral, anticancer* 
antifungal, or other coxxrpound or composition useful for the treatment of a disease. 

"Therapeutic agent regime" is the course of action or use of a therapeutic agent 
35 or combination of merapeutic agents in treating a patient including^ for example, at 
least one of dosage, schedule of administration, choice and/or combination of 
therapeutic agents. 
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'Therapy" is the trealmrat of any disease or abnormality; medical treatment of a 
disease by specified TW>qrw > such as drugs, treatments, or any procedure to ameliorate & 
disease, "^erapy resistance, " as used herein, pertairjeto tiiecar»intyofresista2KX, 
sensitivity, susceptibility, or effectiveness of a therapy against a disease. 
5 'Trough levd" Qoug, is the lowest concentration of a drug in a patient sample 

upon the course of therapeutic agent regimen, 

"Therapy effectiveness" means having the ability to delay progression of at 
l eas t one disease and/or to alleviate at least one disease. 

10 In one emlxxrhnent, one objective of the development of population phannacokmetic 
models for each therapeutic agent is to be able estimate individual phannacoMaetic 
parameters during therapy using one or more plasma concentrations measured at any 
time point after therapy intake and having information on the dosage regimen and the 
time afto the bist drug intake. ' 

15 Previous research has attempted to navigate effective dosages of therapeutic agents to 
challenge rapidly changing ctiologic agents- While the broad approach of population 
pharmacokinetics (usually defined as the change in time of the conceatralim or nature 
of therapeutic agent(s) in groups of patients having similar characteristics) is a 
technique of long standing (see TJsl Ladder*, J. CBn. Pharmacol 28:1059-1062 

20 (1988)), it fills to take into fKTOTiot a large amount of inter-, and even intra-, patient 
variability, ultimately contributing to therapy failure. In diseases such as AIDS, 
therapy failure l eads to the development (and possible eventual dissemination into the 
population) of therapy resistant virus strains. Since neither constant therapy monitoring 
nor completely population-based pharmacokinetic approaches solves all of these 

25 inherent shortcomings, a system and method for optimizing therapy is needed 

Hie problem can be best outlined on thebasis ofan exampla Si^ 
HIV-infected pgfrfryit* receive the same antiretroviral therapy in the same dose three 
times daily. The average pfa?")a concentration-time profile of the therapy in the patient 
population may look as shown in Figure 1 (bold line). However, due to the inter- 

30 individual variability of pharmacokinetic processes (absorption, distribution, 

elimination)* individual plasma concentration-time profiles may subs tantially differ 
fiom the typical profile, as exemplified by the dotted line. A plot of all individual 
plasma concentration-tune profile may cover a range marked by the vertical bars. 
While mdrvidual patient MECs (dashed horizontal lm^ 

3 5 with individual plasma concentraJtran-time profiles or the average plasma 

concenrration-tirne profile, they may cover an area as broad as the grey area Figure 1. 
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As a consequence if the thorny oanaeatntiio& ma patients drops bdow their MEC 
resistance may result 

Because of the importance of maintaining an effective concentration to avoid the 
5 development of disease resistance, and the need to consider an individual patitaifs 

resistance to known therapies in the calculation of optimal dosag e for thqt patient, there 
exists a need in the art for a therapeutic procedure to aid doctors when optimizing 
treatment of these diseases. 

In one embodiment, the present invention, avoids previously known pitfalls in the art 
10 by combining techniques and reiterating obtained data into a model, in order to refine 
the overall model by reducing errors and to generate an or^imized phannacokinetic 
model. This optimized pharmacokinetic model is able to correspond to an individual 
pari cot at a given time, and may be adjusted to correspond to future points in time. 

in another embodiment, the methods of the invention may be adaptable to single 
15 therapies as well as to combinations of therapy regimens and may provide a model with 
inputs for actual individual patient data as well as overall population data from patients 
or individuals (such as from clinical trials), in order to assess for at least one therapy 
whether plasma levels remain above the MEC throughout therapy on a patient by 
patient basis. 

20 In one onbodimeat of the inventi on, the models of the present invention may change 
with time according to the patients' disease progression, new or discontinued drug 
therapy or sensitivity, etc. Systans and methods consistent with the mvention may 
combine at least one bioanajytical method for measuring actual drug concentration in a 
patient at a givoa time, resistance data of the individual pari cut's etiological n Q. pn ^_ and a 

25 first population rAawni^ km^tic model which may foctodtf any relevant fxjvxcriates. In 
one embodiment, the first j b?H 1 fl fM3 fti f T n «+te model may include as much individual 
patient data relevant to treatment as possible to generate dosage($) for all drag(8) which 
will maintain a desired trough level, .above the MEC, for each drug in each patient 
throughout the dosage regimen, whether or not such drugs are currently administered to 

30 the patient 



The systems and methods of the invention may also, for example, include a database 
Corresponding to the data collected and generated from combined first pharmacokinetic 
models and/or from combined optimized pharmacokinetic models. This database may 
35 include a relational geaotypeffienalype database. In a further emrxjdiinaat, a neural 
network or computerized platform may also be provided that learns from the patterns in 
the data collected and generated. 
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In one embodiment of the invention a bioanalytical mefhod is used in the optimization 
of the pharmacokinetic modeL A bioanalytical metbod that may be used in the present 
invention includes, but is not limited to, liquid chromatography with mass spectrometry 
5 (LOMS), An example thereof is a liquid chromatography and mass spectrometry 
assay system currently available from Virco (Mechelen, Belgium), VIRCO 
plasmagram™- Moreover, any other bioanalytical method or methods that provide(s) a 
quantitative measurement of an actual concentration of at least one administered drug 
may be used in the practice of the invention, though bioanalytical methods which 

10 provide a quantitative measurement of all known drags in one or two procedures in a 
short amount of time would provide greater efficiency than methods which require 
longer times and/or more steps. One of skill in the art would realize that in addition to 
the above method, other bioanalytical methods might also be used, such as straight or 
reverse phase liquid chromatography (high pressure or ambient pressure), gas 

15 chromatography, FPLG p rep arativ e chromatography, gel chromatography, ion 

exchange chromatography, etc*, and by detecting with any known detection method, 
suck as fluorescence, UV-vis, IR, NMR, two dimensional multi-wavelength detection, 
etc 

20 For example, in one embodiment, a bioanalytical method may be combined with at 
least one first pharmacokinetic model in order to optimize individual therapy. 
Comparison of the theoretical concentration from the first pharmcokxnetic model and 
the actual concentration is a measure of die accuracy of the first phannacakmetic 
modeL The difference between the theoretical concentration and the actual 

25 concentration may then be minimized by c hangi n g at least one parame te r in the model 
Examples of such parameters include any individual patient data, volume of 
distribution, flhgnrptfon Tflte flftP F frF*. Hft™™ 1 ** ™ ***** tsonstantj etc In one 
embodiment, when the difference is minimized, the pharmacokinetic model is 
optimized for that pati ent at that tune, 

30 

Tn one embodiment, the op timiz ed pharmacokinetic model may be wed in winch at 
least one of three different types of variation and their associated errors are checked 
and minimized: (1 ) -mtra-mdiyi^i iO Vffri fl+K 1 " 1 , where ft single patient's parameters mpey 
change over time (tins includes measurement and g»Hi p l i" g errors); (2) inter-individual 
35 variation, where an individual patient's parameters differ from the calculation baaed on 
previous research and experience; and (3) residual errors, where the theoretically 
predicted drug concentration differs from the actual measured blood drug concentration 
errors. The invention may, for example, address all three sources of error by iterative 
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use of tfaephaTTnnriQknieticmodeL The methods of the invention may also be 
encornpassed uj 4 database, a neural network relating to the database* and/or by the 
combined phannacokinetic model generated fiom previously collected and iterated 
patient data (including data from previously conducted clinical studies), in one method 
S of the invention, a neural network is used to obtain resistance data from genotypic data. 
In another embodiment, a neural network is used to refine the final pharmacokinetic 
model in order to minimize the difference between the theoretical drug concentration 
and the actual cx)ncentraticttL 

1 0 The methods of the invention may also provide, for example, the optimization of 
therapy for a disease such as cancer and/cr retroviral infections (including HIV 
infections in humans or other mammals). The invention also provides a method of 
designing a therapy for a patient, and a method of prescribing a therapy for a patient, 
including making recommendations for drugs and/or combinations of drugs not yet 

15 proscribed for that patient 

Population p ftj™»*^ lcmeric rrj cy ftfl ing 

Population pharmacokinetic modeling is well known in the art Karlsson MO, Shemer 
LB-, J Phannacddnet Biopbarm 1993,21:735-750; Mandema JW, VerottaD, Shemer 

20 LB., J Pbanfcacddnet Biophann 1 £92,20:51 1-528; Any population pharmacokinetic 
modd known in the ait is applicable in the methods of the invention. In one 
embodiment of the present invention, the concentration data obtained by trioana)ysi& of 
human blood samples drawn from a patient is used to develop a population 
: phann a col on eric model. Other information which may be used in such a model 

25 includes, but is not limited to infonnation regarding dosage regimen (dose, dosing 
frequency, therapy formulation, time of administration etc.), the associated sampling 
time, co-medication, and patient-apecific information. 

In one embodiment, a structural pharmacokinetic model may be used in the -p^h^*? of 
30 the invention, which describes the concentratiorrtime course of a therapy. Hie data will 
determine which structural pharmac okin etic model may be used to mathematically 
describe the observed conccntrarion-time courses, 

A population pharmacokinetic model maty describe both the pharmacokinetics of a 
35 therapy in an 'average' patient and the variability of certain parameter values in the 
patient population. 

In population pharmacokinetic modeling, the observed therapy conreaoto&ons in the 
blood may be subject to three types of variability. These are the mter-individual and 
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inter-occasion variability in the pharmacokinetic parameters, mi a residual intra- 
patient, variability. The residual variability originates £rom error in the bio-analysis, 
misspodficatiog of the time After the last drag intake, modd missnenifications etcetera. 
The inter-occaaian variability of model parameters can originate from several causes, 
5 such as variability in hepatic metabolism, increased heart rate, increased water retention 
etcetera- Inte-individual variability ofphannacokindic parameters also originates fiora 
several sources, like the individual's composition of metabolizing enzymes, protein 
composition of the blood, and many others. 

10 A population pharmacokinetic model may comprise co variatcs that explain variability 
of the parameter values. For example, the bodyweigfit of the patient may be predictive 
for a certain pharmacokinetic parameter value for that patient in one embodiment, the 
developed model may be used to predicted pharmacokinetic parameter values of an 
individual patient using Bayesian methods. The obtained parameter values may, for 

15 example, be used to predict the concentration-time course of the drug in that particular 
patient 

In most population based model, fee principal variables are dependent on the model 
used. For example, if a one-<*nnpartmeni model is used, one of the variables may 

20 concern the distribution volume* Since it is difficult to sample a whole patient 

population 24 hours a day, a limited set of sample data is usually available for each 
patient However, the higher the number of patients flic better the estimate of the 
different pharmaco3dnetic variables. In one embodiment, using a given a set of data 
which accurately characterizes the population of interest, the population 

25 pharmacokinetic variables can be readily estimated using software like NONMEM. hi 
another embodiment, the data should consist of a sufficient number of patients to 
riharnrtfriyfl th* p'Karma^Vinf^r. imrUbiKty which grists in die population. Has may 
include deciding which patients to include to cover die natural variability. For example, 
one may include patients in a broad range of weight, age, renal function. 

30 

The NONMEM m* yH for example, provides a quantitative view of the influence of 
various factors including pathological and physiological factors on the 
pharmaccldmetics of the d*ng i «l tti^ population pharniacokmetic par a me ters. Briefly, 
fractional data from individual patients eg. a drug level, may be used to derive 
3 5 population pharmacokinetic parameters which may then be used to derive individual 
patient parameters (viaBayesian approach) again using fiactional data (ag, age, -•-) 
from different individual patients. The patient specific parameters may then be used to 
calculate, for an individual patient, the through concentration or to recalculate the drug 
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dpflaflft to be ttd^faigfr™* to a parinnt fa on? mrKndrme nt, this approach may be used 
to optimize the therapy regimen of an individual patient For example, one may apply a 
Bayesian single c ompartm ent model, . 

5 ^^torvOuotigqt 

As described above, the influence of resistance tearing or therapeutic drug monitoring 
on clinical outcome has often been evaluated separately. Integration of the two areas 
has led to the introduction of a new parameter, tie "inhibitory quotient* (IQ) as a 
potential predictor of clinical outcome. 

10 

The IQ refers to a meagre of the exposure to a therapy in an individual patient (tor 
example; the mipfaxn m concentration, or Qrocgh) divided by die viral susceptibility 
to that therapy in the same patient (for example, IC50 or "fold change" of IC50 as 
compared to wild-type viruses measured in a pheriotypic assay). Other measures of 

1 5 therapy exposure include, but are not limited to, area under the curve, clearance, and 
distribution volume. In one embodiment, the resistance may be determined via a 
VIRTUAIJHENOTYPE® and the virtual IC50 can be used, e.g., IQ may be referred 
to as virtual inhibitory quotient (VIQ). As used here, IQ includes V1Q. Theoretically, 
die IQ orVIQ may be a better measure of resistance because viral resistance is relative 

20 to therapy exposure 

In one embodiment, by relating individual drug exposure to the level of resistance of 
the etiological agent in that same individual, a more accurate prediction of response to 
that drug may be achieved For exanqrie, patients may have ade qu a t e drug levels bat 
25 their etiological agent is moderately resistant, thus they would fail therapy despite good 
drug exposure. The IQ provides additional information over either test alone 
(pheno type or therapy level) and may, &r example, provide clinicians a guide for 
dosage adjustment to achieve Hie desired drug level that can overcome a resistant 
etiological agent 

30 

The normaH2fi»i j^jjKi^ T™rtf¥* 

In one embodiment, the normalized inhibitory quotient (NIQ) is a tool to predict 
clinical outcome using the concept of the inhibitory quotient like the inhibitory 
quotient (IQ), the normalized inhibitory quotient (NIQ) is a ratio of a measure of 
35 therapy exposure and a measure of viral susceptibility to that therapy. However, the 
NIQ corrects for protein binding and may be expressed as follows: 

IQptn ~ IQ of an individual patient determined, for example, by using the actual trough 
concentration and the individual susceptibility of an etiological agent to a therapy: 
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trough concaatralion in the individual paiieat 

IQptn= • ««*- l > 

. fold change of IC50 of the virus in the patient as compared with wild-type virus 

5 

The value of IQpta may then related to the reference intribitoiy quotient (IQref), which 
is the IQ of a patient population. For example, IQref is the mean trough coDcentrabon 
of the therapy as known in the population of patients treated with mis therapy or the 
threshold value for the trough concentration divided by the mean fold change of the 
10 ICSO of a wild-type vims (unity per definition) or the cutoff value of the fold change 
for the normal susceptibility range: 

mean trough concentration in the population or threshold concentration 

IQref- — 

15 fold change of IC50 of wild-type 

Finally, the normalized inhibitory quotient is calculated as follows: 
IQptn 

20 N1Q = <«™ 3 ) 

IQref 

The NIQ may also be multiplied by 100. 

The IQ value provides a direct measure of the success of a patient's therapy In 
25 general, the higher me IQ value, the greater the probability that the therapy is effective. 
Accordingly, the higher the NIQ, the higher the probability that therapy will be 
successful In one embodiment, the NIQ should be around 100%. For example, if me 
NIQ exceeds 100%, the therapy does not need to be changed. While, if the NIQ is 
below 100%, »»m»iid fc* iwvianL either bv increasing the therapy dosage, or by 



30 shifting to a different therapy or a condtination therapy. In one ernbodhnent, the IQ and 
NIQ provide the physician with a single value indicative of the therapy effectiveness. 

Thus, once the IQ is known for at least one tiien^eutin ag«mt, fiw exannile, me 
effectiveness of the at least one therapeutic agent is known and at least cme therapeutic 
35 regime may he optirnized by based on the effectiveness of the at least one therapeutic 
agent Abo, a dosage regime may be adjusted and/or determined, for example, since 
once the IQ is known f» at least ore 

dosage of the at least one therapeutic agent is, for example, known. 

40 
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Adjustment o f Ae dosage irymien for an indtvidual 
In one embodiment, a Bayeaian model may bo used to optimize a population 
phannacokiaetic inodeL The concept of Bayeaian parameter estimation in the field of 
therapeutic drug monitoring is known in the art and may be useful in cijtaimstanccs 
5 where drug concentrations axe measured during relatively complicated dosage 
regimens, or where only a few concentration measurements arc acceptable. The 
Baycsian method allows an estimation of a patient's pharmacokinetic parameters, so 
that therapeutic regimens can be adjusted to achieve specific target concentrations. For 
this purpose, pre-existing information cm population characteristics (means and 
10 variances) of pharmacokinetic parameters is used in conjunction with the (limited) . 
concentration-time data of an individual patient Hie principle of Bayeaian es tima t ion 
is depicted in flow diagram below. 

Individual data £om patients 
(several sources) 

i 

Popolationpharmacokinetic 
modeling 

Population pharmacokinetic 
model 



Data o f individual patient 


— ► 


Bayeaian parameter 




estimation 








* 




Individual phmnacoldnetic 








Target concentration 




Individual dosage regimen 


r — ^ 
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In oidcr to obtain effective treatment, the exposure to the drug (trough concentration, 
AUC, other) sh o uld be higher ft**™ a certain leveL This levd is determined by the 
nature of the etiological population. An indication of the necessary levd maybe * 
5 obtained alter isolating at least one etiological agent and determining the resistance of 
at least one etiological to at least one therapeutic agen^Antrrirogram®, 
Vn-tualPhenotype™, other). 

Generally t phenotypic assays directly measure the ability of a virus to grow in the 
10 presence of each drag of interest, where there may be at least one therapy. One 
technique currently in use, is the ANTIVIROGRAM® (Virco NV, Mechelen, 
Belgium), which is an assay for high-throughput analysis of clinical samples that 
permitts simultaneous detection of HIV phenotypic resistance to both RT and PI (K_ 
Hertogs et aL, Antimicrobial Agents and Chemofcerapy, 42(2): 269-279 (199S), the 
15 disclosure of which is hereby incorporated by reference)- In one embodiment, a 
resistance assay allows an initial estimation of MECs of all known therapies in each 
patient 

The systems and methods of the invention may be implemented through any suitable 
20 combination of hardware, software and/or firmware. Various system components and 
analytical tools, such as neural networks or artificial intelligence, can be utilized to 
farther optimize a drug therapy for the trealment of a disease. In addition, consistent 
with the principles of the invention, a database can be generated through a rombinatioa 
ofbio analytical, population pharmacokinetic, and resistance te s ting methods to provide 
25 individualized therapy regimens that can be administered by physicians and the like. 

The invention may be embodied, for example, as a method, a data processing system, a 
computer program product, a business memod, or any c ombin a ti on thereof Although 
the invention may be practiced without a computer or software-based platform, using a 

30 computer or software-based platfbnn may be desirable, given the complexity of the 
combination and the volume of data ofbioanalytical, population phartn acokin eti c, and 
resistance data obtaining methods. Accordingly, the principles of the invention may be 
implemented as a hardware embodiment, a software embodiment, cs any combination 
thereof and maybe stored many computer usable storage medium, ie%, hard disks, CD- 

35 ROMs, optical storage devices, magnetic storage devices, etc. 

The invention, in one aspect, ia described wife reference to the a c c omp an y ing 
drawings, which include flowchart illustrations of methods and computer program 
products, as well as system or apparatus diagrams. Each block of the flowchart 
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illustrations), or combination of blocks in the flowchart illustration^), cm be 
jmplem miBd by computer pro gram instructions. These computer program instractions 
may be provided to a special purpose computer, a general purpose computer (i.e., a 
computer not dedicated to the methods of the invention alone), or any other data 
5 processing apparatus, to produce a machine such that the ipstnirrions, which execute 
via the processor of die computer or data processing apparatus, create means for 
implementing the functions specified in the flowchart block or blocks. 

Figure 2 provides an exemplify flowchart for optimizing drug therapy, hi one 
1 0 embodiment^ the various steps and operations of Figure 2 may be performed by the 

therapy optimization system 40 in the system environment ofPigure 3 to treat a patient 
diagnosed, for example, with HIV. As indicated above, one of ordinary skill in the art 
will recognize that the features of the ex emplary embodiments can be implemented for 
the treatment of other diseases, such as cancer, other malignancies, or any disease state 
1 5 mediated by a rapidly mutating etiological agent 

As illustrated in Figure 2, in one embodiment the process starts with the gathering or 
collection of patient data (step 100). Patient data may be collected by a physician, a 
doctor or another entity (including clinicians, health care providers, etc). The patient 
20 data may also include the patient's actual drug concentration for one drug, or as many 
dra^ as the patient is taking at that time, and resistance data that is determined from a 
patient sample taken at, or dose -to, that time. In one embodiment, all of the gathered 
patient data may be stored in a database, such as local database 46 of therapy 
op tim iza ti on system 40 (see Figure 3). 

25 

As part of computing an opt imized drug therapy, clinical data is also gathered (step 
110). As part of fttis step, therapy optimization system 40 may include data from 
previous studies (from the same laboratory, and/or from available literature studies) 
and/or from previous patients with the identified disease or condition. The ctinicaJ 

30 data, which, £br example, may be accessed from local database 46 and/or public 

databases) 52, may include data from previous visits from the same patient as apart of 
the clinical data set. The rf»Jjfr?f *r1 data may also include data concerning known inter-* 
drug interactions, such as additional sensitivity or synergy, and known drug 
resistance/piiencty^ Clearly, the order of data coflection is 

35 trrelevar^ arKltbeo TUs patient data 

and clinical data, and any, known conflations between, for exinnple, drugs and 
therapies, maybe included in a first pharmacokinetic modeL 
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This pharmacokinetic modd niay be used to generate a theoretical dnig concentration 
(step 120). The modd may also b c used to determine a theoretical concailratian of any 
drug currently taken by fee patient One embodiment of flic present invention uses a 
single compartment Bayesian modd. 

5 

As illustrated in Figure 2, the theoretical dfug concentration, obtained fenn the 
pharmacokinetic model, and the actual drug concentration, measured fiom the patient 
sample, may then be compared to determine what difference (if any) exists between fee 
theoretical and actual ixrocentrations (step 130). This difference is a measure of model 
10 accuracy. Based on this comparison, a determination is made as to whether the 
difference is minimized (step 140). 

ti 

If the difference is not minimized (step 140% tbm at least one parameter may be 
adjusted in the model (step 150). In one embodiment, the adjustments to the 

15 parameters are m ad e so that the difference between the measured and theoretical 
concentrations is minimized* After adjusting the parametas, the modd calculation 
may be run again to determine a new theoretical concentration (step 120), and the 
process is iterated again (steps 130-150) until the difference is determined to be 
minimized (step 140; Yes). In one embodiment, after minimization, the model maybe 

20 deemed to be a final pharmacokinetic model, optimized for that particular patient at 
that point in time. 

An optimal drug dosage may also, for example, be calculated for that patient at that 
point in time. In one embodiment, the particular patienfs drug concentration should 
25 remain above the minimum effective concentration (step 160). In order to accomplish 
this, the optimized pharmacokinetic model may be used to provide an optimal dosage, 
by changing the actual dose and/or its frequency. 

The information may then be transmitted back to the physjeaan, including 
30 recommendations for dosage increases, decreases, or drug changes. Based on fee 
modcl > which contains information fr o m other clinical studies, and on the patients 
resistance profile, an initial estimation may also be made» optimized for that particular 
patient, as to appropriate dosages for other drugs not yet prescribed to that patient. 

35 Figure 3 is an exemplary system e nv ix tnun entin which the features and methods of the 
invention may be implemented (for example* the methods as shown in Figure 2), As 
illustrated in Figure 3, a communication ch a*™* 30 is provided for facilitating the 
t rans fer of data between various system components and entities. These components 
and entities mchifte one or more physicians 12A-12N who interact with or treat patients 
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(not shown), one ox more laboratories 24A-24N, a therapy qpthruzatian system 40, and 
one or more pubEc databases 52, 

Communication c faaond 30 may be implemented through any single or combination of 
5 channels that allow communication between different people, computers, or locations. 
Hie communication channel may be any system that allows communication between 
the different entities illustrated in Figure 3- 

Each of the physicians 12A-12N collects data for each patient or patients, wherein 
1 0 such data is submitted for analysis by therapy optimization system 40 and/or 

laboratories 24A-24N. The patient data gHfhered by the physicians 12A-12N inchides 
any relevant medical data for thai patient and the patients etiological agent and disease 
or condition, or at least as mudi information as is available^ As illustrated in Figure 3, 
this data can be transfered from each of the physicians 12A-12N to each entity through 
15 communication channel 30. 

During a patient visit, at least one patient sample may be taken by the doctor or other 
entity. The patient sample is salt to one of the laboratories 24A-24N to dctemxine data 
for that patient sample. The patient sarnplemay be obtained at any time, either 

20 concurrently or at a different time as a patient visit and may be provided by a doctor, 
or may be obtained by another professional at a different time and forwarded to die 
appropriate site, such as a laboratory. The data from the sample incl u des the 
concentration of any drugs currently being taken by the patient for the disease or 
condition, and the resistance characteristics of the etiological agent This data may be 

25 obtained from a single- sample or from multiple samples, depending on the etiological 
agent and the drug being taken. The drug concentration and resistance data may be 
provided as part of the patient data to die therapy optimization system 40. 
i 

Therapy optimization system 40 may be implemented through any suitable 
30 combination of hardware, software and/or firmware. For example, therapy 

optimization system 40 ma y be implemented through the use of a personal computer, a 
woriring s tati on, a server or any other computing platform. Software or programmed 
instructions may also be provided for controlling the operations of the computing 
platform, consistent with the principles of the invention. As illustrated in Figure % 
35 therapy optimization system 40 may also include a local database 46 for storing patient 
data. Local database 46 may al so store c li n i r * >1 data or such clinical data may be 
accessed from one or more public datab a&es 52 by therapy optimization system 40. 
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Consistent with the methods of the present invention, merapy optimization system 40 is 
configured to optimize and provide a drop therapy for patients treated by physidans 
12A-12N. As further described below, the optimization of the drug therapy may be 
achieved through a combination ofbioanalyticaL population pnannnrolrine t i c, and 
resistance testing methods to provide individualized therapy regimens that can be 
administered to the patient by a physician. The opthnizid drag merapyniay besentby 
system 40 to physicians 12A-12N in numerous formats (eg., written report, electronic 
file, graphical display, etc.) and may be provided to physicians on fee basis or as a fiee 
or ancillary service. 



In order to demonstrate embodiment* of the invention , an example is presented which 
describes the cptinnzation of treatment of HTV. For example, the methods of the 
invention may be useful in regard to both Pi's (protease mhibimre) and.NNR.TT (ooo- 
nucleosidc reverse tiaiiscriptase inbibitora). OneofslcinmtheortwiBrecogm^mat 
15 the present invention can also be used in connection with the treatment of other 
diseases, and mat various modifications can be made (snch as the use of a neural 
network) m order to optimize therapy for individual patients. 

Tram ple 1 : development ftf * P^" atinn ha3ed pharmacokinetic method 
20 General ouflinfe of an example, methodology 

The data obtained from the quantitative analytical method, i.e., the actual drug 

ovulatory concentration levels, were inputted into a m at hem ati c al model. This model 

was men used to predii* the ormcarttan^ 

prediction, usirtg the model, to^ 
25 sanmliiig, and other asarnm^ 

introduced and/or adjusted to dc« the gap fbimd between me premctrf 

value found through the quantitative analytical model. Validation of the model occurs 

by approximating these variables as cloaely as possible. 

30 A classical population pbaonacolrinetic model may be used to predict an individual 
plasma concentration of a drug using a set of mathematical equations. One em bodime nt 
oftherrresemhrventio^ 
to this modd, at the steady sta^ 
can be expressed as follows: 



35 



f( ? D tV ^IlRjt^ «»f *-■/'»). , expf^ tf ) 
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where Cy is me plasma canMnhntion measured in a patient j at time U- Dj is a 
™»^»n,m.v> dose HAniiristered with an mterdosc interval % P/ symbolizes a set of 
individual phannacokinetic parameters: F), and (volume of distribution, 
5 absorption rate constant and elimination rate constant, respectively)- The latter is equal, 
by definition, to meTatio CL/Fj, where CLj is an individual value of drag clearance. F, 
is a fraction of the dose absorbed after oral administration. It is usually assumed to be 
equal to one, and tons, estimates of CL, and are actually the ratios of clearance and 
volume of distribution to the fiaction absorbed. ^ is a random error reflecting a 
10 residual part of the variability in measured concentration not explained by the model. It 
can often be approximated by the assay error. 

k^CLf and Vjm&y be estimated in each subject and for each drag used to treat this 
patient This is a difficutt task plasma samples to be 

1 5 drawn from a patient It may be substantially simplified if we know the distribution of 
parameters in die patient population! 

20 CLjJ'CEfaJl+ndj 

where ^ Katul CM without subscript^ 

patient population, Often one or more typical pharmacokinetic parameters of a 
particular drag are dependent on patient covariates tike body weight or body surface 
25 area, age, gender, etc. Individual covariates for the patient./ are symbolised by % ^ 
and Gajfor F and C2^ respectively, ift* ffc, and Tfcu are residual variabilities in 
individual V and CL, respectively, wbicii remain uncaqplained after mcrading 
covariate effects in die model. 

30 The population model of a therapy may be known if typical values of each parameter 
are known (in the form of equations that relates them to significant covariates, if any) 
such as residual variabilities in parameters in the pat ent population and a residual 
random error in tile concentration. 

35 Developinp population nharmacokraetic models 

Papulation models for most of the 15 anriretro viral drugs currently used in the 
treatment of HTV^nfected patients have been established: Zidovudine, Lamrvudine, 
pidanoeine, ZaMtabme, Stavudine, Abacsvir, Nevirapme, DtdavTrdioe, Efirvirenz, 
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Saquinavir, Ritonavir, indinavir; Ndfinavrr, Lopduavir, Amprtsnavir. The models fcr 
the remaining drugs may be established u?mg plasma concentrations measured in 
Tjflticrtts dunnc fa oa t m ent (therapeutic drug monitoring data). Also, the population 
model s taken from the literature may be verified/validated using the methods of the 
inventkm. lite population phannaookin etic pr o gram NONMEM based cm the approach 
known as non-linear mixed effect modelling may be used for such modelling- • 

Since several antiretroviral drugs will be administered to each patient, and he or she 
may also receive other drugs like antibiotics, antimycotics, etc* an essential aspect of 
the population model development is searching for drug-drug interactions. If the 
infraction exists it may be included in a model as a covariate. 



Individual prediction using Bavesjan feedback 

Therapeutic drug monitoring usually assumes taking one or two plasma samples per 
IS patient which is not sufficient to find individual estimates of pharmacokinetic 

parameters of the drugs of interest The Bayes appro arii uses both individual plasma 
concentration measurements and popul ation typical values of pharmacokinetic 
parameter together with the variability parameters. Bayesian estimates of individual 
parameters for the patient j\ Pbj, are those which minimise the following objective 
20 function: 



OBJ, 



where the summation is performed over &U concentration measurements and model 
25 p^tiivftBra. <7* is the vari^C 0 of Twngfcpd enttf m the measured co n ce n tra t ion, of a drag* 
12 2 is a set of variances corresponding to interindrvidoal variability in parameters (17)- 6 
is a set of all covariates affecting pharmacokinetic parameters. 

Having Bayesian estimates of individual parameters it is easy to calculate the trough 
30 level by applying the pharmacokinetic model equation again. Moreover, we may also 
accomplish the inverse task: the calculation of the dose magnitude which w01 manriain 
a desired trough leveL Una can be achieved by solving numerically the following 
equation with respect to Djl 



35 C^-faj.D^O 
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where Cfeta^y is infect a nrimnnmi cgxtivttooiiccnfrBtiQfl as estimated by 
Antrviiogmm®. The dose correction accoitfin^ to Bayesian individual predictions is the 
essence of die Bayesian feedback method of therapy individualisaiiosL 

5 The standard Bayesian feedback method described above sometimes results in too high 
maintenance dose exceeding die maximum tolerable dose for a given drug. To avoid 
toxicity one can minimise the difference -/(p #J ,Z>,,*) upon condition < 

Aw- The intodose interval t can also be shortened to avoid toxicity, however, more 
frequent dosing usually leads to poorer compliance This constrained feedback may 
1 0 substantially reduce the risk of drug-related side efforts, however, it may also decrease 
the therapeutic outcome^ 

Example 2: Calculation of Inhibitory quotient 
' Two studies demonstrate the use of thelQor the WQ for the protease mklntois 

15 lopixmvirandind^ ha one study in 56 multiple Pi-experienced, 

NNRTI-*iaive patients treated with lopinavir plus efkviienz and 2 NRTLs, a correlation 
was found between the lopinavir IQ and the % of patients with viral load below 400 * 
copiesfmL at week 24. Hie % of patients with viral load below 400 copies/mL atweek 
24 was 70, 80, and 100% if the lopinavir IQ was <4, 4-15, or > 15, respectively. When 

20 using the lopinavir trough concentration alone, no correlation with virologic outcome 
was found. 

In another study, a VIQ for indinavir > 2 was the strongest predictor of virologic 
response over 48 weeks in patients who failed an ixidmavir^ootaining regimen. 10. In 

25 this study, patients fiul^ 

to a ritonavirfndixfevir 400/400 mg bid regimen* with continuation of the NRITs during 
the first 3 weeks. Thereafter, NRHs were allowed to be switched. Virologic response 
was defined as having a decline of 0-5 log viral load from baseline, or a viral load 
below 50 copiesAnL. The IQ was a better predictor of response than number of 

30 mutations and virtual phenotype fold resistance. 

Tablet. 

Summary of available data on the correlations between IQ or VIQ and clinical 
outcome. 



Drag 


patients 


definition of response 


cut-off 


cori'ttcliun 
factor* 


ra£ 


lopinavir 


52* 


% of patients below 400 copiesAnL 
at week 24 


IQ>15 


0.07 


9 
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prntknts 


definition of response 


cat-off 


correction 
factor* 


re£ 


indinavir 


24** 


% of patients bdow SO copies/mL 
at week 48, or with at least 0-5 log 
drop Sum basdine 


VIQ>2 


0.053 


10 



# correctif«i factor (ICjo °f wild-type virus in the presence of 50% human serum) 
that is multiplied with the fold-change in susceptibility (compared to wild type 
virus) of the viral strain isolated from the patient 
5 * IQ available for 52 out of 56 patients 

** VIQ available far 24 out of 37 patients 



Example 3: Normalized IQ 

This example demonstrates how the normalized IQ may provide information regarding 
10 efficacy of a therapeutic agent The first 2 columns of Table 2 represent the titmgh 
concentration and fold change of the virus for saquinavir. The next 2 columns 
represent what a phannacokinetic modd or resistance testing would advise based on 
these tests alone. The last 4 oofanms represent what a normalized [Q would advise 

based on 4 different scenarios for calculating normalized IQ: 

15 ' 

Method 1: threshold trough / mean fold change wild-type 
Method 2: threshold trough / cut-off fold change 
Method 3: mean trough in population / mean mid change wild-type 
Method 4: mean trough in population / cut-off fold change 

20 

Table 2 



Trough in 


Fold 


Phann 


Virologic 


Me&od J 


Method 


Method 


Mellio 


np/mL 


change 


Model 


advice 




2 


3 


4 


500 


2.0 


^faintain 


Sensitive 


125% 


313% 


50% 


125% 


200 


1.0 




Sensitive 


100% 


250% 


40% 


100% 


500 


5.0 


Maintain 


Resistant 


50% 


125% 


20% 


50% 


1000 • 


5.0 


Maintain 


Resistant 


100% 


250% 


40% 


100% 


100 


OS 


Increase 


Sensitive 


100% 


250% 


40% 


100% 


200 


5.0 


X^Camtam 


Resistant 


20% 


50% 


8% 


20% 


50 


1.0 


Increase 


Sensitive 


25% 


63% 


10% 


25% 


200 


2.5 


J^ifemteifl 


Sensitive 


40% 


100% 


16% 


40% 
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In this example, an IQof around 100% provided evidence that the therapy was 
effective. Furthermore, a decline in IQ indicated that the therapy was becoming less 
effective, while an increase is IQ may indi cats that the drag level is raisin g to toxic 
levels. 

5 

E ttamp jg. 4: (>ptjpu^g Cancer Therapy 

One step for the optimization of cancer therapy is obtaining an actual drug 
concentration. This may be obtained from any patient material which is am enable to 
the bioanalytical method chosen. Examples of samples may be solid or liquid, and may 

10 be excreted and coUecie^ or may be removed firnn the patient Further examples of 
suitable samples include (but are not limited to) biopsies from bone, muscle* organ, or 
skin tissue; fecal, saKva, blood, or tear samples; tumor samples fiura breast colon, 
uterine, prostate, or other malignancies 

The resistance data is also collected, whereto the rmnurnrm effective 

15 concentration (MEC) for at least one drug is determined This data may come from a 
phenotypic assay, i.e,, from testing of any patient derived product that enables the 
determination of MEC of at least one drug against the cancer. 

Alternatively, or additionally, the resistance data may be obtained from 
genotypic data- One method is to s e q uence the genotype, using any one of the methods 

20 well known in the art, and to derive resistance data from a genotype/phenotype 
relational database^ The sequencing can be accomplished on all or a part of the 
genotype, and may focus on a particular oncogene or segment of the genome of 
particular interest, Le., on a known tumor suppressor gene such as p53» 

The method continues siroj^dy to that A first pharmacokinetic 

25 model is used to generate a theoretical drug concentration, which is then compared to 
the actual drug cxrac^tration for that drug in that patient at the specified time. The 
difference between the two concentrations is then minimized by adjusting at least one , 
parameter in the first phaimacoldnetic model. Once the difference is minimized, then 
the pharmacokinetic model is deemed optimized for that patient This optimized model 

30 is then used in combination with the MEC in order to produce an optimized therapy via 
dosage recommendations. 

Example 5: NIQ as a ftcftctor . of ViitjQg ic O qtean& 

HIV resistance testing provides information to clinicians regarding the susceptibility of 
35 a patient's HTV-1 to a drug compared to susceptibility of a reference strain. Although 
this has been shown to predict outcome in salvage therapy, it is unable to provide an 
estimate of whether the patient's drug levels are high enough to inhibit a wild-type Of 
partially resistant st rain Given the wide variability in protease inhibitor concentrations 
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and the common use of phannacokineric boosting to achieve higher concentrations, a 
measure that focoipotatM bom ^ 

susceptibility of themfecting vhm. may be useful in predicting antiviral outcome. This 
example demonstrates the correlation of NIQ with clinical outcome in treatment- 



Methods 

Inclusion criteria included; adults (> 1 8 yrs) infected wife HTV-1 as determined by 
ELISAwith confinnatoiy Western blot; a plasma viral burden of > 500 RNA copies/ml 

10 by bDNA method at a screening visit while receiving a protease inhibitor as a part of 
combination therapy for me preceding 20 weeks with no protease mhibitor drag change 
or dose interruption for > 3 days in the most recent 1 2 weeks; a negative serum or urine 
pregnancy test on the day of enrollment; and a history of no intolerance of ritonavir or 
nelfinavir. Patients were excluded for pregnancy or lactation, prior exposure to 

15 abacavir, amprenavir or efevirenz, concomitant therapy at entry with corticosteroids in 
other man replacement doses, chemotherapy, or mvesngational agents, active, mitteated 
opportmristic infection or other major illnesses, malabsorption or other gastrointestinal 
dysfunction which might interfere with drug absorption or render the patient unable to 
take oral medication, a history of serious rash (erythema multiforme or Stevens- 

20 Johnson syndrome) caused by nevirapine or ddavirdine, or concomitant therapy with 
other drugs that would affect cytochrome P450 metabolism 

Patients were enrolled into three parallel treatment groups mat included abacavir 300 
mg bid, ampremtvir 1200 mg bid, and efeviraxz 600 rag daily with either low dose 
25 ritonavir at 200 mg BID, high dose ritonavir at 500 mg bid, or nelfinavir 1250 mg bid. 

Genotyping (VircoGEN IP*, VIRGO) and VIRTUAL PHENOTYPE™ were 
performed on hasebhe samples. Vhal load data were collected at baseline (mean of two 
pre-therapy samples) and at week 24. Serial phaiinacokinetic samples were coDected 
30 over 12 hours after week 3 for ritonavir-boosted regimens and after week 2 for 
nelfinavir-boosted regimens. 

Ainprenavrr conomtrations in plasma were determined by a vahdated LC-MS/MS 
method. 

The normalised iidtibitory quotient (NIQ) waa deteamhed as: 

NIQ= IP patient 

IQ reference 
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Where the IQ in an individual patient (IQpatient) was m\niiaif>A as ratio of Hie patient's 
tough concentration (Cmin) to the snsce&Mity of tie patient' b vires to the drag, 
expressed as fold change compared to wild type vims (Virtual Phenotypc). ThelQpt 
5 was then, related to the reference inhibitory quotient (IQref), in which fee mean 

population trough concentration of the drug from the product label was divided by the 
cut-off value of the fold change for susceptible viruses. 

For amprenavir, nelfinavir, and ritonavir, the concentration 12 hours after dosing was 
10 used as the Cmux For each drug, relationships between viral load change at week 24 
and the Cmin, fold-change in resistance, andNTQ were fit to a sigmoids! maximum 
effect mode]. 

Seventeen patients were available for analysis with pharmacokinetic data, resistance 
testing, and virologic outcome data at 24 weeks. There were nine patients in the 
nelrrnavir group, four in the low dose ritonavir group* and four in the high dose 
ritonavir group, 

FhannaooJkmctica. 

As shown in Figure 4, the axrrprenavir (AFV) NIQ correlated with outcome at 24 weeks 
(p<0.05). A decrease in viral load to < 400 copies/ml at week 24 was seen in 7/8 
patients achieving NIQ > 3.0 for AFV and 1/9 patients with NIQ < 3.0 (p = 0,003), 
Cmin or phenotype alone were less predictive of outcome than the NIQ for AFV. 
Medians and ranges for Cmin, phenotype and NIQs are shown in Table 3. NIQ values 
for APV were a median (range) of 2.8 (03-41.1). 



Table 3. Individual panapgfoffi 



35 



Drag 




VrrtualPhenotypc 


NIQ 




(ngfal) 


(fold-change) 




Amprenavrx 


1266 


4.0 


2.8 




(264-3453) 


(0.6-8.9) 


(03^1.1) 



All data reported as median (range) 



Example 6: Optimizing 
A. Overview 



WO 02/23186 



-31 



15 



20 



25 
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The invention of opthmzrng a therapy as practiced hor«in for HIV involved a 
series of iterative steps by which individual patient data and overall population data are 
combined and interrelated, which produced the most accurate dosage levels for an 
individual patient Ultimatdy, the inventive process is also able to predict accurate 
5 individual dosage levels for drugs not yet administered to that patient 

lie first step was patient intake, where a complete medical history and 
description were obtained from each patient During this intake step, a patient blood 
samp le (either plasma or whole blood) was obtained, wherein the blood sample 
contained the HIV virus. The intake interview also obtained patient specific data 

1 0 The blood sample or plasma was divided into ah quote for resistance typing of 

the HIV virus and quantitative analysis of the drag levels present in the blood. The 
virus was inactivated prior to being typed. While the viral resistance typing may be 
accomplished by phenotypic or genotypic analysis, or a combination thereof; one 
example is as follows: 

15 B. Viral Resistance Typing: 

Generally, phenotypic assays directly measure the ability of a vims to grow in 
the presence of each drug of interest, where there may b e one drug, or many drugs. 
One technique currently in use, Virco "s ANTTVIROGRAM® (Virco NV, Machete*, 
Belgium), was the first recombinant virus assay for high-thrwghput analysis of clinical 

20 samples that pennitted simultaneous detection of KV-1 phenotypic resistance to both 
RT and PI (BL Hertogs et at, Antimicrobial Agents and Chemotherapy, 42(2): 269-279 
(1998), the entire disclosure of winch is hereby incorporated by reference). Briefly, the 
assay utilized PCR amplification of a fragment of the viral genome obtained from a 
patient's blood sample. The amplified fragments and a pro viral clone lacking the 

25 fragment were de<arTjporated into CD4+ 4 KTT4 cells. Successful combination of the 
provirufi and the amplified fragment within the cells resulted in a recombinant virus 
with a complete HIV-1 genome. This recombinant virus was then grown in cell culture 
to obtain a recombinant viral stock of known concentration. Susceptibility testing of 
the recombinant viral stock in the presence of various antiviral agents and a detection 

30 system based on green fluorescent protein determined which agents inhibit replication 
of the recornbinani virus as of the time that the sample was taken. 

This assay all owed an initial estimati on ofMECs of all known auhretro viral 
drugs m each patient This began the rrrocess which enabl ed (i) selection of most 
effective cambmajion of drugs to be used in the patient and <n) therapy optimization 

35 using a combination of tire patients drug resistance, bibaualysis of drug levels, and 
pharmacokinetic modeling. 
C. Bioanalysis of drug levels 
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Either concurrently or subsequently, another aliquot of the sample or plasma 
was analyzed for levels of all drug? currently administered. One assay method fbr the 
quantitative determination of plasma levels of all autixetro viral drugs in a sample has 
been developed and validated and is detailed below. This procedure is advantageous 
5 because Hie sample volume required was as little as 100 microliters, and the complete 
analytical run could be completed in 15 minutes or less. 

This study validated methods for the quantitative analysis of ritonavir (RTV), 
indinavir (ID V), saquinavir (SQV), nelfinavir (NFV)> nc^rirapine (NVP), ddavirdine 
(DLV), DMP-266 (DMPX ampreoavir (AMV), abacavir (ABV), zidovudine (AZI), 

10 didanosine (DDI), stavudine (D4T), zalcrtabine (DDC) and lamivudine (3TC) in human 
plasma with LCrMS/MS. This embodiment illustrates a single quantitative analysis 
method, though any quantitative analytical method known in the art may he used. This 
quantitative analysis determined the levels of those substances in plasma samples of 
HTV-patients as a part of therapeutic drug monitoring. 

15 Experimental Methods: 

Hie foil owing data and conditi ons validated the detection process for one 
bioanalytical process which may be used according to the invention. The process was 
based on LC/MS, and its accuracy was confirmed for aU relevant storage conditions, 
quality control parameters, etc as follows; 

20 HPLC and mass spechu ni e ti ic conditions 

For practical reasons, two different LC-MS/MS methods were applied for 
quantification of the test substances. The test substances wee divided in two groups 
(group 1 and group 2) dependent on the suitability of analytical methods. For each 
group of test substances a method was validated. ' 

25 Group 1 HPLC and MS-conditions (RTV, IDV, SQV, NFV, NVP, DLV, DMP, and 
AMV): 

The LC-MS/MS conditions for the analysis of the test substances in human plasma for 
Group I were as follows. 

The HPLC Column and Guard Column were both SYMMETRY C18 50 mm x 2.1 mm; 
30 dp=3.5 urn (Waters} (except the guard column was 1 0mm), arid the LC method was run 
at ambient temperature with a flow rate of 0-3 ml/mitL The mobile phase was a 
gradient of Solvent A: 10/90 methanol/Milli-Q; 2-5 mM ammonium acetate absolute 
and Solvent B: 90/10 methanol/NfilK-Q; 2*5 mM ammonium acetate absolute, 



Timefnrinl 


%A 


%B 


% water 


% methanol 


0 




37.5 


60 


40 


0.5 


«,5 . 


37.5 


60 


40 
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Tnnefminl 


%A 


%B 


% water 


% methanol 


0.51 


313 


68.5 


35 


65 


2 


31i , 


68.5 


35 


65 


4 


0 


100 


10 


90 


6 


0 


100 


10 


90 


6.1 I 


62.5 


37.5 


60 


40 


10 


62J 


37.5 . 


60 


40 



10 



Detection: API 300 mass spectrometer (PE-Sciex, Toronto, Canada) 
Interface: Turbo Ionspiay; positive mode; Temp 400°C; flow 5000 ml/W 
Musses monitored Period 1: 
NVP: 266.8 ^ 226-2, 
Dwell time: 1 350 nt% Pause time: 50 ms 
P<*k>d 2: 

DLV: 457.3 2203, SQV: 6713 570.1,IDV: 614.5 421.0, NFV: 568-5 
330.0, RTV: 721 .5 ^ 295.8, DMP: 31&2 243.9, AMV: 506*4 245 J, 
Dwell time: 150 jns, Pause time: 50 ms 
Split ratio no split Injection volume: 3|U 



15 



20 



Gxoxtp 2 HPLC and MS^ondrtions (ABV, AZT, DDI, D4T, DDC and 3TC): 
TheLC-MS/MS conditions for the analysis of the group 2 test substances in hrrman 
plasma samples were as follows. The HPLC column was SYMMETRY CJ8 1 50 mm 
x 3.0 mm; dp^5 jxm, and the guaid column was SYMMETRY C18 20 mm x 33 mm; 
dp=S nm, both from Waters Corporation, Milfaid, MA, USA), TheLC wasnmat 
ambient temperature, with a flow rate of 04 ml/min. Hie mobile phase was a gradient 
of Solvent C: Milli-Q waler with 25 mM ammomnm acetate, and Solvent D: 100 
methanol with 2.5 mM ammonium acetate, according to the table to follows: 



Tfminl 


%C 


%D 


%water 
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70 


30 


70 


30 
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60 


40 


60 


40 
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40 


60 


40 


60 


8,1 


70 


30 


70 


30 


12 


70 


30 


70 


30 



Detection: API 300 mass spectrometer (PE-Sciex, Toronto, Canada) - 
Interface: Positive Turbo ionspray, Temp 350°C; flow: 4000 ml/min 
Masses monitored Period 1: 
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DDI: 229^ -> 11 1.9, D4T: 225,2 -> 127«1, DDC: 211.8 111.9, 3TC: 237.0 -> 

all dwell time: 2S0 ms, pause time: SO ms 
Masses monitored Period 2: 
5 AZT: 268.4 127.1, ABV: 287.4 191.0, 
all dwell time: 600 ms 
pause time: 50 ms 
Splitratio approximately 1 2 (flow to the MS about 130 pl/min) 
Injection volume: 50 |il 

10 

STOOL AND STANDARD SOLUTIONS 

Stock solutions of all test substances of stoop 1 at 1000 jig/ml (weight 
corrected for purity) were prepared by dissolving an exact amount of approxim ately 1 
mg of test substances in methanol. Methanol was added to obtain exact concentrations 
15 oflOOOug/ml. 

Stock solutions of aU test substances of group 2 at 1000 ug/ml (corrected 
weight for purity) were prepared by dissolving an exact amount of approximately 1 mg 
of test substances in methanol Methanol was added to obtain exact concentrations of 
1000 ug/ml. 

20 For each test substance, two stock solutions were prepared, one for the 

preparation of calibration standards (stock solutions 1) and one for the preparation of 
Quality Control samples (stock solutions 2). The stock and standard solutions (working 
solutions, K-references and spike solutions) were stored in the freezer at about -20°C. 

25 CALESB ATTQN STA^Aflpg 

Working solutions containing afl test substances per group were prepared by 
dilution of the corresrxmding stock solutions 1 . The working solutions were used to 
prepare plasma calibration standards by adding 1 vo tume of working so hition to 1 0 
volumes of plasma. The concentrations of the test substances in the working solutions 

30 that were used for validation are outlined in Table 1. 



Ta ble 4 Test substance concentrations (gg/jgj) jg plasma calibration standards 



Name 


Calibration standard reference number 

1 U | 3 U 


Is' 


l» 




Group 1 
















NVP 


100 


200 


500 


1000 


2000 


5000 


10000 


DLV 


82.6 


206 


619 


1858 


4128 


8256 


16512 


IDV 


100 


200 


500 


1000 


2000 


5000 


10000 
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DMP 


100 


200 


500 


1000 


2000 • 


5000 


10000 


RTV 


100 


250 


750 


1000 


2500 


7500 


15000 


SQV 


43.7 


87.4 


218 


655 


1747 


4368 


8735 


NFV 


50 


100 


250 


750 


ZUOU 




1 ruwi 
1UUUU 


AMY 


25 


75 


225 


675 


2000 


5000 


10000 


Group 2 
















DDC 


0J 


0.6 


1.25 


2.5 


5 


7.5 


10 


DDI 


25 


50 


100 


250 


500 


750 


1000 


3TC 


25 


75 


200 


500 


1000 


2500 


5000 


D4T 


25 


50 


100 


250 


500 


750 


1000 


AZT 


25 


50 


100 


250 U 


500 


750 


1000. 


ABV 


50 


100 


250 


750 


2000 


5000 


10000 



The plasma calibration standards were processed according to the work-up 
procedure as outlined above. 

Spike solutions for group 1 were used to prepare pools of plasma, quality control 
samples for group 1 by adding 1 volume of spiking solution to 10 volumes of plasma. 
The spike solutions for group 2 were used to prepare pools of plasma quality control 
samples for group 2 by adding 1 volume of spiking solution to 20 volumes of plasma. 
The concentrations of the quality control samples for each test substance are given in 
TableS. 



10 TableS 



Name 


| Low 


Mid |h&i 


Group 1 


NVP 


120 


1000 


9000 


DLV 


99.1 


16S1 


14861 


IDV 


120 


1000 


9000 


DMP 


120 


1000 


9000 


RTV 


120 


1500 


14000 


SQV 


52.4 


874 


7862 


NFV 


60 


1000 


9000 


AMY 


30 


1000 


9000 


Group 2 


DDC 


036 


2 


9 


DDI 


30 


200 


900 1 


3TC 


30 


500 


4500 ! 
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D4T 


30 


200 


900 


ACT 


30 


200 


900 


ABV 


60 


1000 


9000 



After preparation, the QC^solutions woe aliquo ted and stored at -20*C until use. 
All QCs were processed according to the woik-up procedure as outlined in fee 
experimental part. 
5 K-REFERKNCES 

K-referenoe solution per group consisted of a mixture of all test substance in 
mobile phase at a concentratiaa level of the middle QC. 
VAJjmXm PRQCBDVRE 

For both groups, three analytical batches were processed. Each batch consisted of 
10 » Duplicate set of calibration standards at each of seven concentrations. One set was 
analyzed at the beginning of the analytical batch, and one was analyzed at the end 
of the analytical batch in order to verify the calibration over the time period for 
sample analysis. The time between the HPLC analysis of the two sets was about 20 
hours, which corresponds with the approximate time required for analysis of the 
IS QCs and about 100 samples. 

- Quality Control samples (QCs) at three levels in triplicate. 

- One plasma blank 

- K-ieferences 

20 SYSTEM PERFORMANCE 

The KL-referenccs were used to monitor the performance of the LC-MS/MS 

system. For mis purpose a K-refereoce solution was injected regularly during each 

analytical batch. The mean peak area and its coefficient of variation were calculated. 

RESPONSE FUNCTION 
25 Peak areas of both sets of calibration standards together were fitted using least 

square linear regression For all test substances the optimal weighing fac t or was 

detenninedL 

sENsnwrry(rxoQ) 

The LLOQ (lower limit of quantification) of the test substances was set at the 
30 concentration of the lowest calibration standard. 

PRECISION AND ACCURACY 

The accuracy was shown to be within the calibration range by the following 

procedure. The regression parameters (slope and intercept) were used to determine the 

sample concentrations and to recalculate the concentrations of the calibration standards 
35 on the regression line (dcteonination of the accuracy within the calibration range). 
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The accuracy was deierjnrined as the perceatage relative error (RE). 

The pCTffr fm,tttrsg of Ife Tr^ftKM^ in teona frfgcenmcy and precision was 
established by analysis of quality control (QQ samples and calculation on the 
calibration curve in plasma. 
5 For each of the three concentration levels, the within-batch and between-batch 
precision and accuracy woe detenrdned from the results of the QC samples. The 
within-batch (n=3) and bctwecn-bateh (n=3 of fee mean within-batch determinations) 
precision wens detennined as the coefficient of variation (CV) of the mean areas; the 
accuracy was determined as the percentage relative error (RE). The within-batch and 

10 between-batch precision and accuracy were also determined in QGs of which only 55 
ixl or 27.5 ui was processed* 

■ The absolute recovery was analyzed by tfye following method. Triplicate QCs 
at each of the levels were worked-Up. Also in triplicate, blank plasma was worked^ 
In the last step of the sample preparation procedure, to 100 pi of the extracted blank 

15 100 pJ of 5 mM KKunomum acetate was added cqntangng the relevant test substances at 
a concentration of two times the theoretical cciicentration in end solution. The absolute 
recovery was calculated by comparison of the peak areas of the QCs with the peak , 
areas of the plasma samples that were spiked after processing fee samples. 

The matrix effect on the I^VMS/MS analysis was determined by analyzing 6 

20 different batches of plasma at the lowest QOlcveL Also, several pools of plasma, 
obtained from MTV-patients were used for this purpose. 

Of each plasma batch, in duplicate blank plasma was processed according to the 
sample preparation procedures. In the last step of the procedure to 100 uj of the 
extracted sample 100 pi of 5 mM anxnrauum acetate was added, rontarmng the 

25 relevant test substances at a concentration of two j times the theoretical concentration in 
end solution tor the lowest QOleveL The areas of the test substances in these samples 
were compared with the areas of the test substances in end solution. 

sPECiFicrrY 

The identity of the groupl test substances; (RTV, IDV, SQV, NFV, NVP^LV, 
30 DMP, and AMV) and the group 2 test substances; (ABV, AZX, DDI, D4T, DDQ and 

3TC) was denxrastrafced by the response under the specific MRM conditions of the 

analyte and by the retention time of the analytew The absence of interference was 

verified by processing blank plasma in each analytical batch. 

STABUJIY ANALYSES 
35 a)Freea»^wstabih^ 

after 1 and 4 freeze/thaw 

thawing for 2 hours at >1 5 Q C 
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b) Stability in human plasma al room temperature and at 4°C in the dark; 
Triplicate QCs at the mid QC level ware processed directly after thawing and after at 
least 24 hours storage. 

c) Stability in human plasma at -20°C; Triplicate QCs at the mid QC level were 
5 processed at several time-points after preparation. At least a 2- week interval was 

monitored. 

d) StabUitymtraraan0aamaat55°Cfor4hours: Implicate QCs at the low, 
middle and high QC level were sent to Viroo on dry ice. 

Samples are handled according to proper biohazard procedures* Le., an authorized 
10 person in a bidhazard lab cabinet unpacked the QC*s» The data on the tubes was 
checked with the data on the accompanying list New tabes were prepared and 
identified. The plasma was thawed and transferred into the new tubes. The caps of the 
tubes were decontaminated with ethanoL The sample wan transferred nkto the mraiTiafiriT' 
and heated at 55°C for 4 hams. The samples were cooled to room temperature and 
15 subsequently stored at -SO °C until they were analyzed Samples were maintained <m 
dry ice daring transfers. 

For reference, an additional set of triplicate QCs at the krw> middle and high QC levels 

was sent on dry ice and stored at about -20 °C. Hereafter, the QC's were returned 

together wi th the heated QCs an dry ice and processed* 
20 e) Stability in end solution at room temperature and at4°C in the dark: 

Triplicate QCs at the mid QC level were processed and analyzed within & hours and 

after at least 78 hours of storage 

f) Stability of stock solutions in solvents at -20°C: UV spectra of all test 

substances were measured on dilutions of the stock solutions in DMSO, methanol, or 
25 Milti^ water at several trme^ stock and the end of 

this validation study. The spectra and the extinction coefficients at the absorption 

maxima were compared. The absorbance A (1%, 1cm) was ca lml nfttd 

While the above method has been quality control validated for a single method, 

i.e., high pressure liquid chromatography combined with mass spectrometry, any 
30 quantitative method which separates, identifies *mH quantifies the drugs of interest may 

be used 

Individual MECs and plasma levels of all drugs so obtained are then utilized in 
and incorporated into a population pharmacokinetic model as described below, making 
possible the forecasting of optimal individual drug dosage via Bayesiari feedback. The 
* 35 optimal dosage is defined as the maintenance dose coupled with the inlerdose interval 
which ensures the trough level of each drug remains above the corresponding MEC, but 
below a minimum toxic level. 
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D. Population PhanDacoktactic Analysis 

The population pharaiacolrinetic models for cadi therapeutic drug or 
a n l ire flro viral compound allowed the estimation of die trough level doting therapy for 
each therapeutic compound, using plasma concentrations measured at any time point 
5 after drug intake. This analysis utilizes both the resistance data and the plasma 
concentrations derived from the initial patient sample, and also incorporates any 
relevant patient data obtained at intake. 

This rocfhodology of utilising the MECs and plasma concentrations in the 
pharmacokinetic model may be best explained by way of example. A large grow of 
1 0 HIV-infected patients receives the same antketxoviral drug in die same dose three times 
daily, yielding an overall typical plasma concentration-time profile of the drag for the 
group as shown by the bold line in Figure 1. Inter-individual variability of 
pharmacokinetic parameters gives individual curves which may substantially differ . 
from the typical profile as indicated by the dotted line. If all individual curves are 

15 plotted, they would cover the range marked by vertical bars. If individual MECs are 
shown on the same graph( where the dashed horizontal line illustrates a single example), 
they will also cover some range, as indicated by shading. Due to cyclic behavior of the 
drug concentration profiles, the drug level in some of the patients may drop below their 
MEC, potentially negatively effecting the therapeutic outcome. 

20 Tlie ANTTVIROGRAM® assay (a high throughput, nxombinant virus assay 

which measures die viral susceptibility of a patient sample to all available antiviral 
drugs) provides an individual MEC, and if a trough plasma level of file drug were 
known (shown as a circle on the plot), the dosage may be recalculated in a simple way 
and then modified to get a trough value which exceeds MEC. However, blood samples 

25 are usually withdrawn at random times, and often sampling times do not coincide with 
the time of taking a drug (a square on the plot), precluding the direct calculation of an 
optimal dosage. However, with a population pharmacokinetic model which includes 
estimates of pharmacokinetic parameters in atypical patient and of the mterindividnal 
variability in these parameters across the patient population, aBayesian approach will 

30 estimate the most probable individual parameter estimates and then the dosages may be 
adjusted so as to maintain the trough level which exceeds MEC for that particular 
patient as described 



35 
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1 . A method of measuring the efficacy of at least one therapeutic agent 
comprising: 

5 detemumng an actual concentration of said at least one therapeutic agent; 

determining a pharmacologic exposure using a population pharmacokinetic 
model for said at least one therapeutic agent; 

determining resistance of an etiologic agent towards said at least one therapeutic 

agent; 

10 determining an inhibitory quotient for said at least one therapeutic agent based 

on said pharmacologic exposure and said resistance, and 

using said inhibitory quotient to determine efficacy of said at least one 
therapeutic agent 

15 2. The method of claim X, wherein said inhibitory quotient is a normalized 
inhibitoiy quotient 

3. The method of claim 1, wherein the pharmacologic exposure is a trough 
concentration 

20 

4. The method of claim 1, wherein the resistance is derived from a pheaiotypic 
dc tsmunHtioD 

5. The method of claim 1, wherein the population pharrnfKyfrmctic model is 

25 chosen from & measured population pharmacokinetic T fl n/ ^ and a predicted population 
pharmacokinetic model 

6* Hie method of claim 1, wherein the resistance is determined from a virtual 
phenotype determination, 

30 

7* The method of claim 1, wherein the optimized pharmacokinetic model 
mfnimiym at least one error selected from inlxa-indmdiial* mter-individnaU and 
residual error. 

35 8. The method of claim 1 , wherein the resistance data is obtained fiom a sample 
chosen from at least one of a plasma sample, a blood sample, a saliva sample, a tumor 
sample, a tissue sample, and a bodily fluid sample* 

9. The method of claim 8, wherein the sample is a viius^ntaining sample. 

40 
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10. The method of claim 9, wherein the vims is a retrovirus. . 

1 1 * The method of claim 10, wherein the retrovirus is Human Immunodeficiency 
Virus (HIV). 

1 2. The method of claim 8, wherein the sample contains malignant cdls. 

13. The method of claim 1 wherein the optimized population pharmacokinetic 
model is optimized using- a Bayesian modeL 

14. The method of claim 1 further comprising determining an optimal dosage for all 
therapies in series of therapies- 

15. The method of claim 1 further comprising entering said inhibitory quotient in a 



15 

16. The method of claim 1, wherein the at least one therapeutic agent is an anti- 
infectious compound. 

20 17, Tlte method of daim 16, wheiem 
retroviral agent. 

18. The method of claim 1, wherein the anti-infections compound is an snti-tumoral 
agent 

25 

19. A method of measuring the efficacy of at least one therapeutic agent 
comprising: 

a) ob taining an actual conccntatioia of at least one therapeutic agent in a 
patient at a given time using a bioanalytical method; 
30 b) ralfnilnting » fluvt^r^l fti>nry" t™ti "r» <""d «t least one therapeutic agent 

in said patient at said tirne using a first population phagnacokmctic model; 

c) oht aining a difference by comparing the theoretical concentration of said at 
least one therapeutic agent with the actual concentration of said atleast one thorapanic 
agent in a patient; 

35 d) Trrinnnfrang the difference hy changing at least one parameter in the first 

population pharmacokinetic model in order to generate an optimized population 
pharmacokinetic model; 

c) obtaining resistance data tram said patient; 

f) detcsinraing tie inhibitory quotient for said at least one therapeutic agent 
40 based on said optimized population pharmacokinetic model and said resistance, and 
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g) using said inhibitory quotient to determine efficacy of taud at least one 
therapeutic agent " 

20, The method of d aim 1 9, wherein said inhibitory quotient is a normalized 
5 inhibitory quotient 

21. The method of claim 19, wherein the inhibitory quotient is determined using a 
trough concentration from the optimized population pharmacokinetic modeL 

10 22. The method of claim 1 9, -wherein the resistance is derived from a phenotypic 
determination. 

ti 

23. The method of claim 19, wherein the population pharmacokinetic mo d el is 
chosen from a measured population pharmacokinetic model and a predicted population 

15 pharmacokinetic model 

24. The method of claim 19, wherein the resistance is determined from a virtual 
phenotype determ ination. 

20 25. The method of claim 19, wherein the optimized pharmacokinetic mode) 
minimizes at least one error selected from intm-individual, inter-individual, and 
residual error. 

26. The method of claim 19, wherein the resistance data is obtained from a sample 
25 chosen from at least one of a plasma sample, a blood sample, a saliva sample, a tumor 

sample, a tissue sample, and a bodily fluid sample. 

27. The method of claim 26, wherein the sample is a virus-containing sample. 
30 28. The method of claim 27, wherein the virus is a retrovirus. 

29. The method of claim 28, wherein the retrovirus is Human Immunodeficiency 
Virus (HIV). 

35 30. The method of claim 26, wherein the sample contains malignant c eflsi 

31, Themetbod of claim 19, wherein the re^^ 
throughput screen. 

40 32. The method of claim 19, wherein the optimized population phtin^ tpnirtT^i n 
model is optimized using a Bayesian approach. 
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33. The method of djrim 19 further cpt np nsrn g dete nnriring an op timal dosage for 
all therapies in series of therapies. 

5 34. The method of claim 19 further c ompr ising entering said inhibitory quotient in a 
computer database. 

35. The method of daim 19 further cennprising the use of fee inhibitory quotient to 
provide advice to a physician wherein said advice is chosen from; choice of at least one 
10' of a therapeutic effectiveness of at least one therapeutic agent and dosage of at least 
one therapeutic agent 

3 6. The method of claim 1 9 , wherein the at least one therapeutic agent is an anti- 
infectious compound. , 

15 

37. The method of claim 36, wherein the anti-infectious compound is an anti- 
retro viral agent _ 

38. The method of claim 36, wherein the anti-infectious compound is an HTita» 
20 tumoral fl gtgjt 

39. A method of optimizing at least one therapeutic agent regime 
determining a pharmacologic exposure using an optimized population 

pharmacokinetic model for at least one therapeutic agent; 
25 determining resistance of an etiologic agent towards said at least one therapeutic 

agent; 

d eternising the inhibitory quotient for said at least one therapeutic agent based 
on said pharmacologic exposure and said resistance, and 

using said inhibitory quotient to optimize said at lea^ one therapeutic agent 

30 regime. 

40. A method for determining a dosage regime for at least one therapeutic agent 
comprising: 

deter min ing a pharmacologic exposure using an optimized population 
35 pharmacokinetic model for at least npg therapeutic agent; 

toennimng resistance of an etiologic agent towards said at least one &erapeutic 

agent; 

detenrnning die inhibitory quotient for said at least one therapeutic agent based 
on said pharmacologic exposure and said resistance, 
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and using said inhibitory quotient to determine a dosage regime for at least one 



41. A method for providing advice to a physician regarding at least One therapeutic 
5 agent for at least one patient comprising; 

determining a pharmacologic exposure using an optimized population 
pharmacokinetic model for said at least one therapeutic agent; 

determining resistance of an etiologic agent towards said at least one therapeutic 

agent; 

10 determining the inhibitory quotient for said at least one thcrapctotic ageat based 

on said pharmacologic exposure and said resistance, 

and using said inhibitory quotient to provide advice to a physician regarding at 
least one therapeutic agent for at least one patient 

15 42. A method for providing a report comprising: 

determining a pharmacologic exposure using an optimised population 
pharmacokinetic model for said at least one therapeutic agent; 

dete rmining resistance of an etiologic agent towards said at least one therapeutic 



20 detennimng the inhibitory quotient for said at least one therapeutic agent based 

on said pharmacologic exposure and said resistance, and 

providing a report regarding comprising at least one entry chosen from the 
inhibitory quotient and in&nnation derived from the inhibitory quotient 

25 43. A report comprising a normalized irihft^ 

44. A computer system comprising at least one database comprising at least one 
inhibitory quotient for at least one patient 

30 45. A method of identifying at least ooe therapeutic agent effective against at least 
one etiological agent comprising: 

determining a pharmacologic exposure using an optimized population 
pharmacokinetic model for said at least one therapeutic agent; 

determining resistance of said etiologic agent towards said at least one 
35 therapeutic agent; 

d e t e rmi ni n g the inhibitory quotient for said at least one. therapeutic agent based 
on said pharmacologic exposure and said resistance, ■ 

and using said inhibitory quotient to identify at least one therapeutic agent 
effective against at least one etiological agent 

40 
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